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1 Introduction 

We will apply Lie group analysis for investigating the system of nonlinear 
equations 

Aipt - gpx - fvz = ipx'^'ipz - Ipz^lpx , (1-1) 

Vt + f'ipz = ^xVz - i^zVx , (1-2) 

H = ipxPz - ipzpx (1-3) 
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used in geophysical fluid dynamics for investigating internal waves in uni- 
formly stratified incompressible fiuids (oceans). In particular, the system 
(1.1)-(1.3) with / = was used in [1] to study two non-unidirectional wave 
beams propagating and interacting in stratified fiuid. An exact solution 
of the same system, again in the case when f — 0, was employed in [2] 
for investigating stability of a singe internal plane wave. Weakly nonlinear 
effects in colliding of internal wave beams were investigated in [3], [4] by 
using Eqs. (1.1)-(1.3) with / = 0. The system (1.1)-(1.3) with / 7^ was 
used in [5] to model weakly nonlinear wave interactions governing the time 
behavior of the oceanic energy spectrum. 

In these equations A is the two-dimensional Laplacian: 

A = Dl + Dl e.g. A^, = ^ + ^^A(AV.), 

and g, f, N are constants. Namely, g is the gravitational acceleration, / is 
the Coriolis parameter. The quantity appears due to the density strati- 
fication of a fiuid and is constant under the linear stratification hypothesis. 

We will show in what follows that the system of equations (1.1)-(1.3) is 
self-adjoint (in the terminology of [6, 7]) and use this remarkable property 
of the system for calculating conservation laws associated with symmetry 
properties of the system (1. !)-(!. 3). 

In some calculations, e.g. in Sections 4.7, 4.5, 4.8 it is convenient to 
write Eqs. (1. !)-(!. 3) by using the Jacobians J{ip,v) = ipx^z ~ 'ipzVx, etc., 
in the following form: 

AV't-yp.-M^ J(V',AV'), (1.4) 
vt + Mz^J(^,v), (1.5) 

Pt-\ ipx^J{ip,p). (1.6) 



2 Self-adjoint ness 
2.1 Preliminaries 

We will use the terminology and the following definitions from [6, 7] (see 
also [8]). 

Let X = {x^, . . . , x"") be n independent variables, and u = {u^, . . . , u"^) 
be m dependent variables. The partial derivatives of with respect to 
are denoted by = {uf}, rt(2) = {'^i^}; ■ ■ ■ with 
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where is the operator of total differentiation with respect to : 

Even though the operators Di are given by formal infinite sums, their action 
Di{f) is well defined for functions /(x, -u, «(!),.. .) depending on a finite 
number of the variables x, u, M(2), . . . The usual summation convention 
on repeated indices a and i is assumed in expressions like Eq. (2.1). 

The variational derivatives (the Euler- Lagrange operator) are defined by 

^ = aS. + E™.-A.^. a=l,...,™, (2.2) 

where the summation over the repeated indices ii . . .ig runs from 1 to n. 

Definition 2.1. The afijomi egwa^ions to nonlinear partial differential equa- 
tions 

Fa{x,u, . . . ,U(s)) = 0, a = l,...,m, (2.3) 
are given by (see also [9]) 

F*{x,u,^,...,U(s),^^(s)) =0, a = l,...,m, (2.4) 

where /i — {/i^, . . . , //"*) are new dependent variables, and F* are defined by 

F^{x,u,^,...,u^s),H(s)) ^ ■ (2.5) 

In the case of linear equations, Definition 2.1 is equivalent to the classical 
definition of the adjoint equation. 
Consider the function 

£ = //^F^(x,w, (2.6) 

involved in (2.5). Eqs. (2.3) and their adjoint equations (2.4) can be ob- 
tained from (2.5) by taking the variational derivatives (2.2) with respect 
to the dependent variables u and the similar variational derivatives with 
respect to the new dependent variables 

a = ^ + E(-l)^^n---A.^r;^, a = l,...,m. (2.7) 



S^a d^a / - ''dfX,^..., 
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Namely: 

5C 



Fa{x,U,...,U(^s)): (2-8) 



— ^ F*{x,u,n,...,u^s),l^{s))- (2.9) 
This circumstance justifies the following definition. 

Definition 2.2. The differential function (2.6) is called a /ormaZ Lagrangian 
for the differential equations (2.3). For the sake of brevity, formal La- 
grangians are also referred to as Lagrangians. 

If the variables u are known, the new variables n are obtained by solving 
Eqs. (2.4) which are, according to (2.5), linear partial differential equations 
(2.4) with respect to /i". Using the existing terminology (see, e.g. [10]), we 
will call /x" nonlocal variables. 

Nonlocal variables can be excluded from physical quantities such as con- 
servation laws if Eqs. (2.3) are self-adjoint ([6]) or, in general, quasi-self- 
adjoint ([11]) in the following sense. 

Definition 2.3. Eqs. (2.3) arc said to be self- adjoint if the system obtained 
from the adjoint equations (2.4) by the substitution /i — u : 

F*{x,u,u,...,U(^s),U{s))^0, Q; = l,...,m, (2.10) 

is equivalent to the original system (2.3), i.e. 

F*{x,u,u, . . .,U(s),u^s)) = . . a = 1, . . . ,m, 

with regular (in general, variable) coefficients 

Definition 2.4. Eqs. (2.3) are said to be quasi-self-adjoint if the system 
of adjoint equations (2.4) becomes equivalent to the original system (2.3) 
upon the substitution 

= h{u) (2.11) 
with a certain function h{u) such that h'{u) ^ 0. 

2.2 Adjoint system to Eqs. (1.1)-(1.3) 

Let us apply the methods from Section 2.1 to Eqs. (1.1)-(1.3). In this case 
the formal Lagrangian (2.6) for Eqs. (1.1)-(1.3) is written 

c = (p [Ai>t - gpx - fvz - i^x'^i'z + ^l^z'^i'x] 

r 1 ^2.12) 

+ ij'[vt + fi^z - i^xVz + i^zVx] + rypt + —i^x- i^xpz + i^zpx\ , 



where /i and r are new dependent variables. The adjoint equations to 
Eqs. (1.1)- (1.3) are obtained by taking the variational derivatives of £, 
namely: 

where (see (2.2); see also Eqs. (3.6)) 
5 
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Taking into account the special form (2.12) of C, we have: 



+ 



— = - -D,---{Dl + Dl) 

5i) dipx dip. 



■ dC ^ dC ^ dC 
DtTT^ + Dx^^^ + Dz 



Dx {ipAijjz + iJiVz - —r + rpz) - (ipAipx + fl^ + l^'Vx + rpx) 



{Dl + Dl) Dticp) + DxicpiPz) - Dzi^iPx) 



7V2 

Vx^'il^z - VzAlpx + l^xVz rx + TxPz - fP'z - l^zVx " rzPx 

9 

- Aift + 2 [ifxzi'xx + '^zzi'xz - ^xxi'xz - ^xzi^zz] , 



SC ^ dC ^ dC ^ dC , ^ 

Dt- Dx^ = -l^t - fJ'xWz + JVz + l^zWx , 

dv OVt OVx OV^ 



5jC ^ dC ^ dC ^ dC 
dp dpt dpx dpz 



-n + g^x - r-x^Jz + r-zipx 



Hence, the adjoint equations (2.13) can be written as follows: 

Aift + — + /At. - VxAipz + VzAiPx - e = 0, (2.14) 

- Pt- fJ'xi^z + f<Pz + fJ'zi^x = 0, (2.15) 

-rt + g^Px - Txipz + Tzipx = 0, (2.16) 
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where 

Q = J{lJ,,v) + J (r, p) + 2 [(pxzi^xx + ^zzlpxz - ^xxlpxz - ^xzlpzz] ■ (2.17) 

2.3 Self-adjointness of Eqs. (1.1)-(1.3) 
Theorem 2.1. Eqs. (1.1)-(1.3) are quasi-self-adjoint. 

Proof. Looking for (2.11) in the form of a general scaling transformation, 
one can readily obtain that after the transformation 

(f^iP, 11^ -V, r^-—p, (2.18) 

the quantity given by Eq. (2.17) vanishes. Therefore the adjoint equa- 
tions (2.14)-(2.16) become identical with Eqs. (1.1)-(1.3) after the sub- 
stitution (2.18). Hence, according to Definition 2.4, Eqs. (1.1)-(1.3) are 
quasi- self- adjoint. Since Eqs. (2.18) are obtained just be simple scaling of 
the equations ip = 'il)^^ = v.r = p required for the self-adjointness, we will 
say that Eqs. (1.1)-(1.3) are self-adjoint. 



3 Conservation laws 

3.1 General discussion of conservation equations 

Along with the individual notation t, z for the the independent variables, 
and p, if) for the dependent variables, we will also use the index notation 
— t, — X, x^ — z and — v, v? — p, — ip, respectively. We will 
write the conservation laws both in the differential form 

Dt{C') + D,{C') + D,{C') = (3.1) 

and the integral form 

j^J J C^dxdz = 0, (3.2) 

where the double integral in taken over the the (x, z) plane M?. The equa- 
tions (3.1) and (3.2) provide a conservation law for Eqs. (1.1)-(1.3) if they 
hold for the solutions of Eqs. (1.1)-(1.3). The vector C = {C^^C^^C^) 
satisfying the conservation equation (3.1) is termed a conserved vector. Its 
component is called the density of the conservation law due to Eq. (3.2). 
The two-dimensional vector (C^, C^) defines the flux of the conservation 
law. 
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The integral form (3.2) of a conservation law follows from the differential 
form (3.1) provided that the solntions of Eqs. (1.1)-(1.3) vanish or rapidly 
decrease at the infinity on M^. Indeed, integrating Eq. (3.1) over an arbitrary 
region f2 C we have: 

^JJ C^dxdz^-J J [D^{C^) + D,{C^)]dxdz. 

According to Green's theorem, the integral on the right-hand side reduces 
to the integral along the boundary dQ oi Q : 

-If [D^{C'^) + D,{C^)\ dxdz = f C^dx - C^dz, 

and hence vanishes as ft expands and becomes the plane R^. 

Remark 3.1. It is manifest from this discussion that one can ignore in 

"divergent type" terms because they do not change the integral in the 
conservation equation (3.2). Specifically if evaluated on the solutions of 
Eqs. (1.1)- (1.3) has the form 

C^^C^ + D^{h'') + D,{h^) (3.3) 

with some functions h'^,h^, then the conservation equation (3.1) can be 
equivalently rewritten in the form (see [12], Paper 1, Section 20.1) 

Dt{C') + D,{C^) + D,{C')^0, 

where 

= C^ + Dt{h^), = C^ + Dt{h^). 
Accordingly, we have 



J J C^dxdz = J J C^dxdz, 



and hence the integral conservation equation (3.2) provided by the conser- 
vation density of the form (3.3) coincides with that provided by the 
density C^. 

In particular, if = the integral in Eq. (3.2) vanishes. This kind 
of conservation laws are trivial from physical point of view. Therefore we 
single out physically useless conservation laws by the following definition. 
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Definition 3.1. The conservation law is said to be trivial if its density 
evaluated on the solutions of Eqs. (1.1)- (1.3) is the divergence, 

C'^D,{h'') + D,{h^). 

The following statement ([13], Section 8.4.1; see also [7]) simplifies cal- 
culations while dealing with conservation equations. 

Lemma 3.1. A function F{v, p,ip,v^,v^, p^, Pz,ipx,il)z,ip^t,'4^zt, ■ ■ ■) is the 
divergence, 

F = D,{C') + D,{C'), (3.4) 
if and only if satisfies the following equations: 

Sv ^ Sp ^ Sip ^ ^ 

Here the variational derivatives act on F as usual (see also Section 2.2): 

Sv dv ^ \dvxj \dvzj 

_dF _^ f9F\_^ fdF\ 
Sp dp \dpxj ^ \dpz) 

5F _dF _ ^ f dF\ _ ^ f dF\ ( \ [ ^F \ 

5^} di) ""ydilj^) KdijjJ \dipxt) ^ \d^zt) 

CoroUeiry 3.1. A function is the density of a conservation law (3.1) if 
and only if the function 

F = DAC^) (3.7) 

^ ' (1-1)-(1.3) ^ ' 

satisfies Eqs. (3.5). Here |(i.i)-(i.3) means that the quantity Dt{C^) is 
evaluated on the solutions of Eqs. (1.1)-(1.3). 

In particular. Lemma 3.1 allows one to single out trivial conservation 
laws as follows. 

Corollary 3.2. The conservation law (3.1) is trivial if and only if its density 
evaluated on the solutions of Eqs. (1.1)-(1.3), i.e. the quantity 

^* ^ ^^l(i.i)-(i.3) (3-^) 



satisfies Eqs. (3.5), 



5Cl 5Cl 5Cl , , 

0, ^ = 0, ^ = 0, (3.9) 



bv ^ bp ' bil) 
on the solutions of Eqs. (1.1)-(1.3) 
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3.2 Variational derivatives of expressions with Jaco- 
bians 

We will use in our calculations the following statement on the behaviour 
of certain expressions with Jacobians under the action of the variational 
derivatives (3.6). 

Proposition 3.1. The following equations hold: 

^^ = 0, ^^^ = 0, (3.10) 
dv dip 



5v 


= 0, 


6[vJ{ilj,v)] 
5ijj 


0, 


(3.11) 


Sp 


= 0, 


6[pJ{ip,p)] _ 


0, 


(3.12) 


SJ{iIj, At/;) 


= 0, 


5[^J(V^,AV^)] 

8il) 


= 0. 


(3.13) 



Proof. Let us verify that the first equation (3.10) holds. We have (see 
(3.6)): 



Sv Sv 



Replacing v hy ip one obtains the second equation (3.10). Let us verify now 
that Eqs. (3.11) are satisfied. We have: 

S[vJ{ij,v)] S[v{ip^v^-ip^v^)] d[v{ip^Vz-ip;,v^)] I \ , r I \ 

dv Sv dv ^ ' ^ ^ 

= J (-0, v) - D^ivil)^) + D^ivil)^) = J{ii,v)- J {ii, v) - vil)^^ + vil^zx = 

and 

S[vJ(aP-v)\ „ / „ / 

— = -D^{vv^) + D^[vv^) = -v^v^ - vv^^ + v^i;^ + vv^^ = 0. 

Replacing v by p one obtains Eqs. (3.12). Eqs. (3.13) are derived likewise 
even though they involve higher-order derivatives. We have: 

(5J(Va At/;) _ Sgj^Ajj^ ~ IpzAll^x) _ Sjll^xjll^xxz + i^zzz) - i^zji^xxx + i^xzz)] 

Sip Sip Sip 

= -L>,(AV'.) + D,{AiP,) - D,{Dl + Dl){iP,) + D,{Dl + Dl){iP,) 
- -D^iAiP,) + D,{Ai/j,) - D,{Ai/j,) + D,{Ai{j,) = 0. 
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Derivation of the second equation (3.13) requires only a simple modification 
of the previous calculations. Namely: 

= ipxAip^ - ipzAipx - ipxAipz - ipAi/j^x + ipzAipx + ipAipxz 
- ^AD.DxiiP^) + ^ADxD.iiP^) = 0. 



3.3 Nonlocal conserved vectors 

It has been demonstrated in [8, 7] that for any operator 

admitted by the system (1.1)- (1.3), the quantities 

dC 



(3.14) 



du 



dC 



dC 

rill" 



DjDk{W-) 



dC 



(3.15) 
i = l,2,3, 



define the components of a conserved vector for Eqs. (1.1)-(1.3) considered 
together with the adjoint equations (2.14)-(2.16). Here 



= ry" - ^J-u" a = 1, 2, 3. 



(3.16) 



The formula (3.15) is written by taking into account that the Lagrangian 
(2.12) involves the derivatives up to third order. Moreover, noting that the 
Lagrangian (2.12) vanishes on the solutions of Eqs. (1.1)-(1.3), we can drop 
the first term in (3.15) and use the conserved vector in the abbreviated form 



du] 
r dC 



(3.17) 



dC 
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For computing the conserved vectors (3.17), the Lagrangian (2.12) con- 
taining the mixed derivatives should be written in the symmetric form 



~ '4^x{^zxx + '4^xzx + '4^xxz + '^'4^zzz) + '4^z{^'4^xxx + "^^xzz + '4^zxz + '^zza;)] 



r 

+ P[Vt + flpz - 'IpxVz + IpzVx] +r Pt^ Ipx- IpxPz + 'IpzPx 

L 9 

(3.18) 

Since the Lagrangian £, and hence the components (3.17) of a conserved 
vector contain the nonlocal variables (fi, p,r, we obtain in this way nonlocal 
conserved vectors. 



3.4 Computation of nonlocal conserved vectors 



The substitution of (3.18) in (3.17) yields: 



dvt 



dpt 
dC 



Di 



dC 



dtp, 



txx 



+ D 



dC 



dtp, 



tzz 



dip, 



txx 



dC 



dip-, 



tzz 



+ Dl{W 



dC 



dip. 



+ Dl{W') 



txx 



or 



(3.19) 



1 
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Furthermore, using the same procedure, we obtain: 



dC 



+ D. 



dC 



+ Dl 



-Dt{W')D., 



dC dC , ^ ^ 
dC 



dC dC 

+ 



dipxtx 
dC 



-Dx{W') Dt 



d^pxxz Oil), 
dC 



dC 



-Dx{W')Dx 



dC 

di)xx'. 



dipxxt 

dC 

di^xx 



-D,{w') a 



dC 

dipxz^ 

dC 



dC 



Olpxzz 



dC dC 

+ 



dipxtx dipxxt 



or 



dil^x 
+ D,D,{W^) 



dC dC 
+ 



di^xxz dip 

X 



= fii/j^ + (ri/j^ -gip) + W^\- Ai/j^ -l^v^ + — r - rp^ (3.20) 

9 

+ Dl{^iP,) + ^ Dli^iP,) + \^xt-\ DxD.i^iPx)] 



+ ^ ipDtDxiW^) - ^ ipiPxD.DxiW^). 



<P1pz 



Likewise we get 

= -W^ (piPx + M - W riPx + [ - A^, + (/ + vx)p + rpx (3.21) 



rl 



1 



DziW-") 77 - D.i^ijjx) + - Dx{^ipz) 
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3.5 Local conserved vectors 



The quantities (3.19)-(3.21) define a nonlocal conserved vector because they 
contain the nonlocal variables </?, /i, r. In consequence, the conservation equa- 
tion (3.1) requires not only the basic equations (1.1)-(1.3), but also the 
adjoint equations (2.14)-(2.16). 

However, we can eliminate the nonlocal variables using the self-adjointness 
of Eqs. (1.1)-(1.3) thus transforming the nonlocal conserved vector into a 
local one. Namely, we substitute in Eqs. (3.19)-(3.21) the expressions (2.18) 
for (p, ii,r : 

(2.18) 
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Then the adjoint equations (2.14)-(2.16) are satisfied for any solutions of 
the basic equations (1.1)-(1.3), and hence the quantities (3.19)-(3.21) satisfy 
the conservation equation (3.1) on all solutions of Eqs. (1.1)-(1.3). 

Let us apply the procedure to C^. We eliminate the nonlocal variables 
in (3.19) by substituting there the expressions (2.18) and write in the 
following form: 



where 



= (^) + Dl (^) , AW^ = Dl{W^)+ Dl {W^). 
We further simplify the expression for by using the identities 



and 



Then we have: 



i,Dl[W^)=D^[iljD.,[W^)]-ij.,D.,[W^), 



C^ = -vW^-l-pW-'-'il^, D,{W'') - D,{W') + - A{,PW') . (3.22) 



Dropping in (3.22) the divergent type term 



-A{^|;W')^D, 
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in accordance with Remark 3.1, we finally obtain: 

C^^-vW^-^pW-i^^D, {W"") - {W"") . (3.23) 

We will not dwell on a similar modification of the expressions (3.20), 
(3.21) for the components and of conserved vectors. We will see 
further in Section 4.5 that they can be found by simpler calculations when 
a density is known. 

4 Utilization of obvious symmetries 

4.1 Introduction 

Eqs. (1.1)-(1.3) do not contain the dependent and independent variables 
explicitly and therefore they are invariant with respect to addition of arbi- 
trary constants to all these variables. It means that Eqs. (1.1)-(1.3) admit 
the one-parameter groups of translations in all variables, 

V = v+tti, p = p+a2, ip = ip+as, t = t+a4, x = x+a^, z = z+Gq, 

with the generators 

Y - ^ Y - ^ Y - ^ Y - ^ Y - ^ 
ov op dip at ox 

One can also find by simple calculations the dilations (scaling 
tions) 

V — av, p — bp, ip = cip, t = at, x — (3x, z — f3z (4.2) 

admitted by Eqs. (1.1)-(1.3). These transformations are defined near the 
identity transformation if the parameters a, . . . , P are positive. The dilations 
of x and z are taken by the same parameter /3 in order to keep invariant the 
Laplacian A. Let us find the parameters a, . . . , (3 from the invariance con- 
dition of Eqs. (1.1)-(1.3). The transformations (4.2) change the derivatives 
involved in Eqs. (1.1)-(1.3) as follows: 

Vt = aavt, = af3v.^, % = a/3vz, 

pt = bapt, Px = bPpx, Pz = bppz, (4.3) 



(4.1) 
transforma- 
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where A is the Laplacian written in the variables x, z. The invariance of 
Eqs. (1.1)- (1.3) under the transformations (4.2) means that the following 
equations are satisfied: 

AV'f - gpx - fvz - i^x^i^z + i^z^i^x = o, 

Vi + - ll^xVz + = 0, 

_ _ 

pt H - i^xPz + i^zPx = 0, 

9 

whenever Eqs.(l.l)-(1.3) hold. Substituting here the expressions (4.3) we 
have: 

AV'f - gpx - fVz - i^x^i^z + i^z'^i^x = cap^Ai/jt - b/3gp^ - c^/3'^{'iIj^A'iIj^ - i/j^Ai/j^), 
vt + fi'z - i'xVz + i'zVx = aavt + c^fip^ - acf3'^ {ip^v^ - ip^v^) , 

Pt H - i^xPz + i^zpx = hapt + 0/3 — - bcp U^p^ - ip^p^) . 

g g 

These equations show that the invariance of Eqs. (1.1)-(1.3) is guaranteed 
by the following six equations for five undetermined parameters a,b,c,a, (3 : 

^bp ^c^l3\ aa^c/S ^ac/S^, ba ^ c/S ^ bc/3^ . (4.4) 

It can be verified by simple calculations that Eqs. (4.4) yield 

a = 1, b = a, c = a^, /3 = - , 

a 

where a is an arbitrary parameter. We substitute these values of the pa- 
rameters in (4.2), denote the positive parameter a by e", drop the tilde 
and conclude that Eqs. (1.1)-(1.3) admit the one-parameter non-uniform 
dilation group 

t^t, X ^ xe", z = ze", V = ve", p = pe", t/j = ^e^" 
with the following generator: 

9 d d d ^ , d ^, 
^7 = 2;— + + t;— + p— + 2^— • 4.5 
ox oz ov op dip 

Wc will consider the operators (4.1)-(4.5) as obvious symmetries of Eqs. 
(1.1)-(1.3) and will compute the local conservation laws provided by these 
symmetries. 
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4.2 Translation of v 

For the operator Xi from (4.1) Eqs. (3.16) yield 

= 1, W"^ = 0, = 0. 
Substituting these expressions in Eq. (3.23) we obtain 

C = -V. 

In this case the equations (3.20) and (3.21) are also simple. They are written 

and upon using Eqs. (2.18) yield: 

C = -vi;,, = vi,^ - fi,. 

Since any conserved vector is defined up to multiplication by an arbitrary 
constant, we change the sign of C^, C^, and obtain the following con- 
served vector: 

C^ = v, C^ = vij,, C' = fij-vij,. (4.6) 

We have: 

A(C^) + D,{C^) + D,{C') ^v, + v.i^, + /V'. - 
Hence, the conservation equation (3.1) coincides with Eq. (1.2). 

4.3 Translation of p 

For the operator X2 from (4.1) Eqs. (3.16) yield 

= 0, W'^ = 1, = 0. 
Substituting these expressions in Eq. (3.23) we obtain 

Furthermore, Eqs. (3.20), (3.21) and Eqs. (2.18) yield: 

Multiplying C^, C^, by —N'^jg^ we arrive at the following conserved 
vector: 

C^=p, C^^— V^ + pV., C^-pi^,. (4.7) 
9 

One can readily verify that the conservation equation (3.1) for the vector 
(4.7) is also satisfied. Namely, it coincides with Eq. (1.3). 
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4.4 Translation of if) 

For the operator from (4.1) Eqs. (3.16) yield 

= 0, = 0, = 1. 

Substituting these expressions in Eq. (3.23) we obtain 

=0. 

Hence, the invariance of Eqs. (1.1)-(1.3) under the translation of ^0 furnishes 
only a trivial conservation law (see Definition 3.1). 

4.5 Derivation of the flux of conserved vectors with 
known densities 

We will show here how to find the components and of the conserved 
vector (4.6) without using Eqs. (3.20), (3.21), provided that we know the 
conserved density = v. 

Let us first verify that — v satisfies Corollary 3.1. In this case 
DtiC'^) = vt, and hence Eq. (3.7) yields 

F = A(C^) ^J{i^,v)-fi^,. (4.8) 

(1.1)-(1.3) 

Using Proposition 3.1 we see that Eqs. (3.5) are satisfied: 

Therefore Corollary 3.1 guarantees that F defined by Eq. (4.8) satisfies Eq. 
(3.4): 

i^.v, - i'zv^ - M = D,{li^) + D,{ii'') (4.10) 

with certain functions if^, . 

In order to find if^, if^, we write 

■^xVz - ji>z = D^{V^^ - f^) - V^:cz, -i^zVx = D^{-V^z) + Vll^zx 

and obtain: 

i^xVz - IpzVx - f'ipz = D^{-V1pz) + D^{vip^ - fljj). 

Thus, = -vijj„ = vi)^ - fijj. Denoting = -H\ = -//^ i.e. 
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and invoking Eq. (4.8), we write Eq. (4.10) in the form 

A(C^) +D,{C') + D,{C^)=0. 

(1.1)-(1.3) 

This is precisely the conservation equation (3.1) for the vector (4.6). Thus, 
we have obtained the components C^, of the conserved vector (4.6) 
without using Eqs. (3.20), (3.21). 

The components C^, of the conserved vector (4.7) can be derived 
hkewise. 



4.6 Translation of x 

For the operator X5 from (4.1) Eqs. (3.16) yield 
Substituting these expressions in Eq. (3.23) we obtain 

Hence, the invariance of Eqs. (1.1)-(1.3) under the translation of x furnishes 
only a trivial conservation law (see Definition 3.1). Similar calculations 
show that the invariance under the translation of z provides also a trivial 
conservation law. 

4.7 Time translation 

For the operator from (4.1) Eqs. (3.16) yield 

Substituting these expressions in Eq. (3.23) we obtain 

^ VVt + ppt + ^x^xt + ^z^zt- (4.11) 

Changing the last two terms of by using the identity 

= L>,(V'V'.0 + ^.(V'V'.t) - ^^A (4.12) 

and dropping the divergent type terms, we rewrite given by Eq. (4.11) 
in the form ^ 

C'^vvt + j^ppt-^ljAijt. (4.13) 
18 



Let us clarify if the conservation law with the density (4.13) is trivial 
or non-trivial. According to Definition 3.1, we have to evaluate the density 
(4.13) on the solutions of Eqs. (1.1)-(1.3). In this case it is convenient to 
use Eqs. (1.1)-(1.3) in the form (1.4)-(1.6) and replace Eq. (3.8) by 

^1 ^1 1 

^* — ^ l(1.4)-(1.6)' 

Then we have 

2 

Cl = {vJ{^, v) + j^ pJiij, p) - ^J{^, A^) } - fD,{vi,) - gDx{pi,) 

and Corollary 3.2 shows that the conservation law is trivial. Indeed, the 
last two terms of C\ have the divergent form. The expression in braces for 
C\ satisfies Eqs. (3.9) according to Proposition 3.1, and hence it also has 
the divergent form. 

Thus, the invariance of Eqs. (1.1)-(1.3) under the time translation fur- 
nishes only a trivial conservation law. 

4.8 Use of the dilation. Conservation of energy 

Consider the generator (4.5) of the dilation group, 

d d d d ^ , d 
ox az ov op aip 

In this case the quantities (3.16) have the form 

— v — xvx — zvz, W'^ = p — xpx — zpz, — 2ilj — xijjx — z'ijjz- (4.14) 

The substitution of (4.14) in (3.23) yields: 

= -v'^ + xvv^ + zvv^ + {-p^ + xppx + zppz) 

-^1 + Xll^^ll^^cx + Zll^xi^xz Xipz'ipxz + Zlpz'ipzz- (4.15) 

We modify (4.15) by using the identities 

1 1 

xvVx + zvVz = -Dx (xf ^) + -Dz {zv^) — v"^, 

Xppx + Zpp^ ^\^^ i^P^) ^\^^ (^^^) ~ 

Xlpxi^xx + Xlpzi^xz ^^D^[x {ipl + - ^ {ipl + Ipl) , 

Z^x^xz + Z^,^,, ^\Dz[z {^l + ^l)] -1(^1+ ^l) ■ 
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Substituting these in (4.15) and dropping the divergent type terms we have: 



where 



Dividing by the inessential coefficient (—2) we finally obtain the following 
conservation law in the integral form (3.2): 



~dt 



dxdz = 0. (4.16) 
Eq. (4.16) represents the conservation of the energy with the density 

(4.17) 



Let us find the components and C"^ of this conservation law written 
in the differential form (3.1). We will use the procedure suggested in Section 
4.5. Let us first verify that E defined by Eq. (4.17) satisfies Corollary 3.1 
for densities of conservation laws. We have 



Dt{E) = 2(vvt + -^ppt + il^x'ipxt + 'ipzipzt ] ■ 



(4.18) 



Since the expression in the brackets in Eq. (4.18) is identical with (4.11) it 
can be rewritten in the form (4.13), and hence satisfies Eqs. (3.5). Corollary 
3.1 guarantees that E is the density of a conservation law. It is manifest 
from Eq. (4.17) that this conservation law is non-trivial. 

According to Corollary 3.1, Dt{E) defined by Eq. (4.18) and evaluated 
on the solutions of Eqs. (1.1)-(1.3) satisfies Eq. (3.4), 



Dt{E) 



(4.19) 



(1.1)-(1.3) 

with certain functions H^, H'^. In order to find H^, H"^, we use Eq. (4.12), 

2(^.^a:t + i^zi^zt) = ^.(2^^a:i) + ^.(2^^.t) - 2ijAi,t, (4.20) 

and write: 



2vvt = ^il^^vv^ - ip^vv^ - 2fvil)^ 



(4.21) 
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7V2 



L>,(p'V'.) - - 2yL',(pV') + 2#p., (4.22) 



= -2g^Pp, - 2fijv, - D, A^,) + {i^^A^l^.) . (4.23) 

Substituting the expressions (4.21), (4.22) and (4.20), (4.23) in the right- 
hand side of Eq. (4.18), we arrive at Eq. (4.19) with 

Thus, denoting = —H^^ = —H^ we arrive at the following differential 
form (3.1) of the conservation of energy for Eqs. (1.1)-(1.3): 

Dt{E) + D,{C^) + D,{C^) = (4.24) 

with the density E given by Eq. (4.17) and the flux given by the equations 

= 2^pV + v'^^l^z + 1^ p'^z - 2^^xt + V'^AV'., 

= 2fvij - ^; V. - 1^ - "^i^i^zt - i^'^Mx- (4.25) 

5 Invariant solutions 

5.1 Invariant solution based on translation and dila- 
tion 

Let us find the invariant solution based on the following two operators: 

d d 

mX^ — kX^ = m— k— (m, k = const.), 

ox oz , , 

(5.1) 

d d d d ^ , d 
ox oz ov op aip 

We first find their invariants J(i, x, z, v, p, ip) by solving the equations 

(mXs - AjXg) J = 0, X7J = 0. (5.2) 
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The characteristic equation kdx + mdz = for the first equation (5.2) 
yields that the operator 171X2 — kX^ has, along with t, v, p, ip, the following 
invariant: 

X = kx + mz. (5.3) 

Therefore we have to find the invariants J{t, A, v, p, ijj) for the operator X7. 
To this end, we write the action of X^ on the variables t, A, v, p, iJj by the 
standard formula 

d d d d 

and obtain 

d d d d 

To solve the equation X-jJ{t, A, p,ip) — for the invariants, we calcu- 
late the first integrals for the characteristic system 

dX dv dp dip 
X V p 2ip 

and see that a basis of invariants for the operators (5.1) is given by 

Accordingly, we assign the invariants V, i?, to be functions of the invariant 
t and arrive at the following general form of the candidates for the invariant 
solutions: 

v = XV{t), p = XR{t), ijj = X'^(p{t), X = kx + mz. (5.5) 

In order to find the functions V (t) , R{t) , (l){t) , we have to substitute the 
expressions (5.5) in Eqs. (1.1)- (1-3). 
We have: 

ijjt = A20'(t), = 2/cA0(t), = 2mA0(t), 
VVt = 2(A;2 + m^)<f>'{t), V^V. = 0, V^V^, = 0, 
i/j^V;, = 2A;mA 0(i) V{t) , t/j^v^ = 2A;mA 0(0 V{t) , 
ip^p^ = 2A;mA (j){t) R{t) , VzPx = 2A;mA 0(t) R{t) . 
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Therefore Eqs. (1.1)- (1.3) yield the following system of first-order linear 
ordinary differential equations: 

2(A;2 + m^)<f)' - gkR - fmV = 0, 
XV' + 2/mA0 = 0, 
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or 



V' = -2fm(l), (5.7) 

i?' = -2-Ar20. (5.8) 
9 

Let us integrate Eqs. (5.6)-(5.8). Differentiating Eq. (5.6) and using 
Eqs. (5.7)-(5.8), we obtain 

(f)" + = 0, (5.9) 

where 

The general solution of Eq. (5.9) is given by 

(j){t) — Cicos{u!t) + C2sm{u!t), Ci,C2 = const. (5.11) 
Substituting (5.11) in Eqs. (5.7)-(5.8) and integrating, we obtain 

2/m 



2k 



R^C^ 

guj 



Ci sin(a;i) — C2 cos(a;i) j , 
Ci sm{ujt) — C2 cos{ut) 



To determine the constants C3 and C4, wc substitute in Eq. (5.6) the above 
expressions for V, R and the expression (5.11) for and obtain 

/mCs + gkCi = 0. 
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Thus, the solution to Eqs. (5.6)- (5.8) has the following form: 

4>{t) = Ci cos{ut) + C2 sm{u!t), 
2fm 



V{t) 

m 



UJ 

2k 

gu 



C2 cos{u!t) — Ci sm{u!t) 
C2 cos{ujt) - Ci s\n{ut) 



+ C3, 

fm 
gk 



(5.12) 
(5.13) 

(5.14) 



Finally, substituting (5.12)-(5.14) in (5.5), we arrive at the following 
solution to the system (1.1)-(1.3): 



guj 



C2 cos{ut) — Ci sin(a;t) 



gk 



2fm 



C2 cos(wt) — Ci sin(a;t) A + C3A, 



-0 = [Ci cos{ujt) + C2 s\n{ut)\ A^, 



(5.15) 

(5.16) 
(5.17) 



where A is given by (5.3), u is defined by Eq. (5.10) and Ci,C2,Cz are 
arbitrary constants. 



5.2 Generalized invariant solution and wave beams 

It is natural to generalize the candidates (5.5) for the invariant solutions 
and look for particular solutions of the system (1.1)-(1.3) in the following 
form of separated variables: 

v^F{\)V{t), p^a{X)R{t), ip ^ p{X)(p{t), X^kx + mz. (5.18) 

The reckoning shows that then the right-hand sides of Eqs. (1.1)-(1.3) 
vanish and Eqs. (1.1)-(1.3) become: 

{k^ + m^)p"{X)(j)\t) - gka '{X)R{t) - fmF\X)V{t) = 0, (5.19) 

F{X)V'{t) + fmp'{X)(f){t) = 0, (5.20) 

kN"^ 

a{X)R'{t) + -J- P\X)(t){t) = 0. (5.21) 

Differentiating Eq. (5.19) with respect to t, using Eqs. (5.20)-(5.21) and 
dividing by we obtain 



0. 
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Assuming that the ratios /3' , a' / a, F' / F are proportional with constant 
coefficients and have one and the same sign, we arrive at an equation of the 
form (5.9). For example, letting 

B" a' F' 

TT = — = — , (5.22) 

we obtain Eq. (5.9). Then, according to (5.11), we can set in (5.18) 0(i) = 
cos{ujt) and 0(i) = sm{ujt), i.e. 

ijj = A{\) cos(wt) and = B{\) sin(wt). (5.23) 

For each function i/j given by (5.23) we determine the functions V{t), R{t) 
using Eqs. (5.20), (5.21), (5.21), then take the linear combinations of the 
resulting functions and arrive at the following form of the "generalized in- 
variant solution" (5.18): 

^ = ^(A) cos{ujt) + B{X) sin(cjt), (5.24) 

v^^— [B'{X) cos{ujt) - A'{X) sm{ujt)] + F(A), (5.25) 

kN^ 

p = [B'(X) cos(a;i) - A'(X) sm(ujt)] + H(X), (5.26) 

guj 

where u is given by Eq. (5.10). 

The reckoning shows that the functions (5.24)-(5.26) with arbitrary 
A(A), B{X) solve Eqs. (1.1)-(1.3) provided that F(A), H{X) satisfy the 
following equation: 

gkH'{X) + fmF'{X) = 0. (5.27) 

One can readily verify that the invariant solution (5.15)-(5.17), which is a 
particular case of (5.24)-(5.26), obeys the condition (5.27). 

5.3 Energy of the generalized invariant solution 

If we substitute in (4.17) the generahzed invariant solution (see Eqs. (5.24)- 
(5.26)) 

V' = A{X) cos{ut) + B{X) sm{u}t), (5.28) 
v^im: [B'{X) cos{ujt) - A'{X) sm{ujt)] , (5.29) 

00 

p = [B'(X) cos(ujt) - A'(X) sm(ujt)] , (5.30) 

gu 
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where 



A — kx + mz, u! — 



we obtain: 

E^{k'' + m')[A'\\) + B'\\)\. 

Invoking that any conserved vector is defined up to multiphcation by an 
arbitrary constant, we divide the above expression for E by (A;^ + m^) and 
obtain the following energy: 

E = A'^{\) + B'^{X). (5.31) 

Since the energy density (5.31) depends only on. \ — kx -\- mz, it is 
constant along the straight line 

kx + mz — const. (5.32) 

Accordingly, the "local energy" (5.31) has one and the same value at points 
(a;o,^o) and {xi^Zi) provided that 

kxQ + mzQ — kxi + mzi. (5.33) 

The energy density (5.31) describes the local behavior of the solutions. 
Therefore it is significant to understand its distribution on the {x, z) plane. 
Suppose that the functions ^(A), B{\) and their derivatives rapidly de- 
crease as 77 — > oo. If we take, as an example, the functions 

where a is a positive constant, then the energy density (5.31) of the wave 
beams has the form 

E 



(1 + A2)2 

Hence, the energy is localized along the straight line (5.33). Therefore we 
can define a wave beam through a point {xq, zq) as the totality of the points 
{xi,zi) satisfying Eq. (5.33), i.e. identify it with the straight line (5.33). 



6 February 2009 



26 



Bibliography 



[1] A. V. Kistovich and Y. D. Chashechkin, "Nonlinear interactions of two- 
dimensional packets of monochromatic internal waves," Izv. Atmos. 
Ocean. Phys., vol. 27, pp. 946-951, 1991. 

[2] P. N. Lombard and J. Riley, "On the breakdown into turbulence prop- 
agating internal waves," Dyn. Atmos. Oceans, vol. 23, pp. 345-355, 
1996. 

[3] A. Tabaei and T. R. Akylas, "Nonlinear internal gravity wave beams," 
J. Fluid. Mech., vol. 482, pp. 141-161, 2003. 

[4] A. Tabaei, T. R. Akylas, and K. G. Lamb, "Nonlinear effects in re- 
flecting and colliding internal wave beams," J. Fluid. Mech., vol. 526, 
pp. 217-243, 2005. 

[5] R. N. Ibragimov, "Latitude-dependent classification of singularities for 
resonant interaction between discrete-mode internal gravity waves," J. 
of Physical Oceanography. To appear. 

[6] N. H. Ibragimov, "Integrating factors, adjoint equations and La- 
grangians," Journal of Mathematical Analysis and Applications, 
vol. 318, No. 2, pp. 742-757, 2006. doi: 10.1016/j.jmaa.2005.11.012. 

[7] N. H. Ibragimov, "A new conservation theorem," Journal of Mathe- 
matical Analysis and Applications, vol. 333, No. 1, pp. 311-328, 2007. 
doi: 10.1016/j.jmaa.2006.10.078. 

[8] N. H. Ibragimov, "The answer to the question put to me by 
L.V. Ovsyannikov 33 years ago," Archives of ALGA, vol. 3, pp. 53- 
80, 2006. 

[9] R. W. Atherton and G. M. Homsy, "On the existence and formulation 
of variational principles for nonlinear differential equations," Stuides in 
Applied Mathematics, vol. 54, No. 1, pp. 31-60, 1975. 



27 



[10] I. S. Akhatov, R. K. Gazizov, and N. H. Ibragimov. "Nonlocal symme- 
tries: Heuristic approach," Itogi Nauki i Tekhniki. Sovremennie prob- 
lemy matematiki. Noveishye dostizhenia, vol. 34, pp. 3-84, 1989. En- 
glish transL, Journal of Soviet Mathematics, 55(1), (1991), pp. 1401- 
1450. 

[11] N. H. Ibragimov, "Quasi-self-adjoint differential equations," Archives 
of ALGA, vol. 4, pp. 55-60, 2007. 

[12] N. H. Ibragimov, Selected Works, III. Karlskrona: ALGA Publications, 
2008. 

[13] N. H. Ibragimov, Elementary Lie group analysis and ordinary differen- 
tial equations. Chichester: John Wiley & Sons, 1999. 



28 



